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c/3 . Abstract 

We present a review of some of the current major challenges in stellar cluster research, includ- 
ing young clusters, globular clusters, and galactic nuclei. Topics considered include: primordial 
■ mass segregation and runaway mergers, expulsion of gas from clusters, the production of stellar 

exotica seen in some clusters (eg blue stragglers and extreme horizontal-branch stars), binary 
populations within clusters, the black-hole population within stellar clusters, the final parsec 
problem, stellar dynamics around a massive black hole, and stellar collisions. The Modest Ques- 
tions posed here are the outcome of discussions which took place at the Modcst-6A workshop 
held in Lund, Sweden, in December, 2005. Modest-6A was organised as part of the activities of 
the Modest Collaboration (see www.manybody.org for further details). 



1 Introduction 

MODEST is an abbreviation for MOdeling DEnse STellar systems, and is a collaboration between 
groups working throughout the world on stellar cluster research, including both theoreticians and 
observers. The Modest-6A workshop was held in Lund, Sweden, in 2005, as part of the continuing 
activities of the Modest Collaboration. A particular task of this workshop was to produce a list 
of challenges in stellar cluster research - The Modest Questions - considering both problems likely 
to be solved in the shorter term (around one year) and those requiring more work (timescales of 
several years). This paper provides a review of what came out of the discussion at Modest-6A. The 
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topics have been grouped into three areas (young clusters, globular clusters, and galactic nuclei) 
although there is naturally some overlap between the three sections. 



2 Young Clusters 



2.1 Primordial mass segregation and runaway mergers 

Dynamical interactions between stars in yo ung clusters leads t o mass segregation, in which the heav- 
ier stars sink towards the cluster centres ( Bonnell fc Davie's! . 19981 : de Griis et al. . 20021 : KroupaL 
20041 : Fleck et al. . 200fJ ). This occurs very rapidly, on a timescale of roughly t mseg ~ ireiax^av/^heavy, 
where t rc \ ax is the two-body relaxation time of the cluster and m av ,7Bheavy are the average stellar 
mass and the mass of the heavy stars; t mseg may approach the crossing time for very young dense 
clusters. Howev er, clusters may form in a mass-segregated state by the competitive accretion pro- 
cess outlined in iBonnell et al,l (|l997n . Whether or not mass segregation is a primordial state of 
star clusters is a crucial question because it affects the timescale for core collapse and it also has 
some bearing on the origin of massive stars. This in turn affects the process of mergers involving 
massive stars in cluster cores. A cluster which is intially mass-segregated is more likely to undergo 
a runaway merger process, since the massive stars are already situtated at the cluster centre. 



In a star cluster where the two-body relaxation time is sufficiently small (t 



relax ~ 



100 Myr) 



the most massive stars can reach the cluster core by dynamical friction and drive the cluster to 
a sta t e of core collapse b efore they explode ()Portegies Zwart &; McMillan! . l2002al : | Giirkan et all 



Freitag et all EoORl ) . During this phase the stellar number density in the cluster core be- 
comes so high (n c ^ 10 8 stars/pc 3 ) that stars may experience direct physical collisions. This 
may lead to a runaway colli sion process in which one star repeatedly collides with other stars 



( Portegies Zwart et al. . fl999h . The growth rate of this object may exceed 10 _3 MQ/yr, and can 
therefore exceed t he stellar mass loss rate, which for the most massive stars is of the same order 
( Vink et alll200ll ). The mass of the single massive object may grow to about 1000-3000 Mq al- 



though the subsequent evolution of this object is unclear; it may for example collapse to form a 
black hole ( Portegies Zwart et al. . 20041 ). Observations may tell us whether or not such interme- 
diate mass black holes form in dense stellar clusters. Go od candidates so far a re the young and 
dense star cluser MGG-11 in the starburst galaxy M82 dMcCrady et all l2003l ) . and the stellar 



conglomerate IRS13E very near the Galactic center (|Maillar d et al J,|2004J). The lat ter object is of 
parti cular interest as its black hole may be of the order of 1500 M to 15000 M ([Schodel et al. . 
200,51 ). The subject of intermediate-mass black holes within globular clusters will be discussed in 
section 3.7. 



For the MODEST questions and tasks on this subject, we propose: 



In the next year: perform N-body simulations of clusters with and without primordial mass 
segregation to determine its effect on the evolution of the clusters, particularly in regard to 
core collapse and runaway mergers. 

In the next three to ten years: Do any mass-segregate d systems exist which are too young to 



m tne next tnree to ten years: uo any mass-scgrcgatcc 
have segregated dynamically ( Gouliermis et all 20041 )? 



In the next ten years: How do winds and instabilities affect the evolution of massive stars 
and what implications does this have for the behaviour of very massive merger products? 
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2.2 Gas expulsion from clusters 



Embedded clusters represent a crucial but poorly- understo od phase in the process by which a 



giant molecular cloud is converted to a population of stars ( Dale et al. . 20051 : Clark et al. . 2005h . 



The observation tha t most embedded clusters do not survive to become open or globular clusters 
( Lada fc Ladall2003h implies either that most embedded clusters become unbound during the star- 
formation proc ess, or that the g iant molecular clouds from which they form were never bound in 
the first place ( Elmegre'enl . 2000h . 

There are essentially three types of embedded systems ( Kroupa fc Boilv . 20021 ): Type / contain 
from a few to 1000 stars but no O stars since they are too rare to be sampled from the IMF. 
Type II have between about 10 3 and 10 5 stars and contain between a few to about one hundred 
O stars. Type III clusters are massive with > 10 5 stars. Gas-expulsion may take a few crossing 
times for type I clusters because the cumulative feedback energy from low-mass stars drives the 
gas out. For type II clusters gas expulsion may be explosive because the O stars provide sufficient 
feedback energy to blow out the gas on a crossing timescale or shorter. In the very massive type 
III systems, feedback from O-stars is likely to be insufficient to remove gas until the detonation of 
the first supernovae because the gas density in these systems is large enough to quench the effects 
of photoionisation and winds. Such massive systems can therefore achieve core collapse while still 
containing substantial quantities of gas and therefore modelling them poses a particularly difficult 
problem as one must allow for both the dynamical N-body and hydro dynamical evolution of the 
system. 

The efficiency with which feedback expels gas determines the star-formation efficiency and also 
the likelihood of the cluster becoming unbound: a cluster unbinds more readily for l ow star for 



mation efficiencies an d /or for gas expulsion timescales shorter than a crossing time ([Lada et al 



198 



Goodwinlll997l 'l. It is therefore important to model this process correctly. Gas expulsion also 



converts a hydrodyna mical problem into an N-body problem, and t hus determines when N-body 
calculations may start ( Gever &: Burkertl . 200ll : Kroupa et al. . 2001 ). At present, stellar feedback 
is poorly understood because it is not possible to treat directed outflows self-consistently, nor is it 
possible to handle the full radiative transport problem in three dimensions. 



We propose the following MODEST questions and tasks: 

• In the next one to three years: creation of a hybrid N-body SPH code, using high-precision 
N-body codes such as NBODY6 or starlab as a starting point, since most existing hybrid 
codes are SPH codes which have had an N-body component grafted on. 

• In the next ten years: What is the efficiency with which O-stars expel gas from typical clusters 
and hence what is the typical efficiency of star-formation? 

• In the next ten years: study the formation of a 10 6 Mq cluster from the collapse of its natal 
molecular cloud until the stage in cluster evolution where mass loss becomes dominated by 
stellar evolution 



2.3 Cluster complexes 

It has been recognized for some time that studying the formation of isolated stars does not paint a 
realistic picture of star formation, since virtually all stars form in clusters and many interact with 
each other during their formation. There is also a considerable body of evidence that cluster forma- 
tion is itself clustered. Interacting galaxies have been known for some time to host vast networks 
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of very ma ssive star clusters whose forma tion is trig gered by mergers or tidal interactions, e.g. the 
An tennae more fc Schweizerl . M 995ah and M82 (jde Oriis et all Emh. However, as pointed out 

by lLarsen &: Richtlerl (|1999h the formation of such clust er associations is not co nfined to merging 



galaxies. Dwarf galaxies (e.g. NGC 1569 and NGC 1705. lO'Connell et al 



spiral galaxies (e.g. NGC 253. IWatson et al.1 dl99fih. M10 1 
large star-forming complexes. Elmegreen fc ElmegreenI 



Bresolin et al 



1994)) and undisturbed 



1996J)) also exhibit very 



J19.8J3) identified regularly-spaced strings 
of giant HII regions or HI clouds in the spiral arms of 22 galaxies, finding that individual concen- 
trations of star-formation had typical size and mass scales of 1-4 kpc and 10 6 -10 7 M respectively, 
considerably larger than the sizes and masses of typical molecular clouds or OB associations. 

Star-formatio n in the Milky Way itself also appears to occur on scales much larger than tradi- 
tionally assumed. Shevchenkol (| 19791 ) identifies more than t en indicators of star formation and lists 
49 star-formation regions with masses of up to ~ 10 6 M . Efremov ( 19781 ) used Galactic Cephied 
variab les to identify 35 complexes of star formation with sizes of ~ 600 pc and iBarkhatova et al 



( 19891 ) used spatial and kinematic data on 69 open clusters to identify eleven star complexes with 
sizes ranging from ~ 10 2 -~ 10 3 pc. Efremo\] (|l979l ) suggests that the sizes of these regions simply 
reflect the original sizes of the complexes of gas and dust from which they formed. He then pro- 
poses that the spatial and temporal distribution of star-formation within the complexes is a result 
of self-propagating star-formation driven by feedback from O-stars. 

The study of star cluster complexes is therefore of great importance, since it will shed light on 
the importance of triggering in star formation on a variety of lengthscales. The evolution of such 
complexes is cl early an important factor in the formation and evolution of individual star clusters 
( Kroupal . 1998|) and may also have a beari n g on t he origins of Ultra-Compact Dwarf galaxies and 

E002al lbT l2005h . 



'faint fuzzies' ( Fellhauer fc Kroupal . 



We propose the following MODEST questions: 

• In the next year: Is star cluster formation itself always clustered, or do some clusters form 
and evolve alone? 

• In the next three years: From an observational perspective, how do cluster complexes evolve 
and what is the role of triggering (on any scale) in their formation? 

• In the next ten years: What bearing, if any, does the evolution of star cluster complexes have 
on the formation of the various dwarf-galaxy populations and the large faint star clusters 
(faint fuzzies)? 



3 Globular Clusters 



3.1 The production and evolution of blue stragglers 

Blue stragglers are main-sequence stars which are more massive than the current turn-off mass. 
They have been seen in globular and open clusters, and in the halo. They are believed to form 
in two ways: either through collisions between lower-mass main-sequence stars (which will be 
relatively frequent in dense stellar clusters (H i lls &: Day, 1976)), or via mass transfe r within binaries 
(Preston fcSnedel Ewl IPiotto et all ESj IPavies et all Eq3; ISandnuistl ^M)- 

The subsequent evolution of merger and mass-transfer products remains an open question. Hy- 
drodynamic calculations suggest that the collision remnants do not develop substantial convective 
regions dur i ng th eir thermal relaxation and therefore are not mixed significantly after the collision 



( Sills et all . 119971 ). If true, this would have a significant impact on the subsequent evolution of the 
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merger product, as the core would receive only a small amount of hydrogen and would, therefore, 
have a relatively short life. However, another key ingredient is the angular momentum contained 
in the merger product. Even a small, non-zero, impact parameter can result in an object with 



substantial angular momentum (jSills et all 120021 ) . This poses a problem for its subsequent evolu 



tion. The total angular momentum can be up to ten times larger than that possessed by low-mass 
pre- main-sequence stars. The collision product is expected to have quite a large radius soon after 
the merger and to then gradually contract back towards the main-sequence. However if no angular 
momentum were removed during the process, the object would reach the breakup velocity long 
before it reaches the main-sequence. Some mechanism must remove most (~99%) of the angular 
momentum. In pre-main-sequence stars angular momentum loss is driven by surface convection 
zones and magnetic w inds, which are not expected to be present in collision products. Recently 
De Marco etakl d2004h inferred from spectroscopic data the existence of circumstellar discs around 
6 stars in a sample of 50 objects located above the main-sequence in 4 globular clusters. The 
presence of magnetic ally locked discs had al ready been suggested as a possible mechanism to lose 
angular momentum ( Leonard fc LivioL 19951 ). The observed discs might not be massive enough to 
explain such an angular momentum sink, but they could be the leftovers of once larger discs. 

Unfortunately very little observational data for BSS rotation rates in gl obular clusters currentl y 



exists. Only a few BSS in globular clusters have measured rotation so far (|Pe Marco et all l2005h . 



Comparison between model predictions an d observat i on can thus be made only on the basis of 
colour-magnitude diagrams. For example, ISills et al.1 (l2005h claim that disc-free models of BSS 
are brighter and bluer than the objects observed so far and therefore imply that some angular 
momentum loss mechanism must be at work. 

We propose the following MODEST tasks: 

• In the next one to three years: spectroscopic observations of BSS will shed some light on 
the origin of BSS and on their subsequent evolution. Rotation rates, surface gravity and 
chemical abundances are fundamental information for probing collision models. Collisional 
BSS are expected to form at the centre of globular clusters. For this reason observing them 
spectroscopically will be quite challenging and will probably require instruments such as 
HST/STIS and ground-based adaptive-optic systems. 

• In the next one to three years: improve models of rotating collisional products, and develop 
more detailed modelling of mass transfer in binaries. Predicted differences in the observables 
can help distinguishing collisional from primordial BSS. 



3.2 Understanding the observed properties of extreme horizontal-branch stars 

Extreme horizontal-branch (EHB) stars have been observed in several globular clusters as a group 
of objects considerably bluer than regular horizontal-branch stars. It is by now widely accepted 
that EHB stars are He burning sta rs that during their evolution have suffered heavy mass loss 
( Iben fc Rood! Il970l : iFaulknerl . Il972h . keeping only a thin envelope (with mass of the order of 0.02 
M Q ). However the actual formation mechanism f or such object is still unclear. Mass loss du ring the 
horizontal branch (HB) phase has been proposed ( Wilson fc Bowenl . ll984l : lYong et al.|. 20001), as we ll 



200ll ). 



as enhanced mass loss rates during the previous red giant branch (RGB) phase (poker et al 
through several mechanisms. One of these in volves binarity , i.e. t i dal interacti o ns within a close 
binar y could enhance the envelope mass loss dMengel et all Il97fil: iHeber et all 120021 : iHanetal 



120021 ), However, preliminary results from spectroscopic data ( Moni Bidin et alll2006h show a lack o f 



ZUUZ). riowever, preliminary results irom spectroscopic data l|lvlom oicun et al.l.lzu un ) sliow a iacK_p 
binaries among EHB stars in NGC 6752, in sharp contrast with the results of lPeterson et ah ( 2002h 
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that concluded that the majority of EHB stars in this clust er are in binary system s. Apparently the 
two catalogues sample different regions in the cluster, the Peterson et alJ (|20f)2l ) one being located 
in more external regions. This could imply a different formation mechanism for EHB in the central 
and in the outer regions as has been observed for BSS in many clusters by looking at t he luminosity 



functions (Ferraro et al 



2004). Recent result s on the very massive clusters, uj Cen ((Piotto et al 



20051 ) and NGC 2808 (jP'Antona et all l2005l ). indicate that the EHB stars in these two clusters 
could be the product of a second generation of stars formed by material enriched in He due to 
the pollution by SN and/or intermediate mass AGB stars. This would also explain the correlatio n 



found between the extension of the EHB and the mass of the clusters ([Recio-Blanco et al.l . 11997 s ). 



We propose the following MODEST question and task: 

• In the next one to three years: could dynamical interactions in dense stellar clusters trigger 
the large envelope mass loss believed to be at the orgin of EHBs? 

• In the next one to three years: understand the binary distribution among EHB stars by 
observing other clusters. 



3.3 Multiple episodes of star formation in some globular clusters? 

Photometric studies of red giants in u Cen have revealed several discrete populations covering 
the metallicity rang e —1.5 < [M/H] < —0.5 with, possibly, an age spread of up to ~ 6 Gyr 



( Sol lima et all |2005j) . Additional observations have revealed the presence of a do uble main se : 



quen ce, with a population of stars lying to the blue of the primary main sequence (jBedin et al 
2004h . Spectroscopic follow-up measurements provided the s urprising result tha t the blue main 



sequence is ~ 0.3 dex more metal-rich than the red population (jPiotto et al. l. l2005l V the most likely 
explanation being that the blue main sequence represents a super He-rich population of cluster 
members. This interpretation is interesting because u Cen also possesses an EHB, for which one 
viable explanation is a population of He-rich stars. 

A number of other globular clusters are also known to possess unusual stellar populations. 
Recent HST observations of NGC 2808 have demonstrated this cluster to also have a population 
of blue main seque nce objects, although the main sequence does not show the clear bifurcation 
present in u Cen ( D'Antona et al. . 20051 ) . Again, this population has been interpreted as He- 
rich - a scenario which may also help explain NGC 2808's EHB. In addition, two Galactic bulge 
clusters (NGC 6388 and 6441) also possess EHBs despite their rather high metallicities ([Fe/H] ~ 
—0.5). Helium enhancement has been invoked to explain the anomalous HBs in these two objects. 
M54, which lies ne ar the centre of the Sagi t tarius dSph galaxy, is suspected to possess a small 
metal licity spread dDa Costa fc Armandroffl. Il995l ). and has recently been shown to possess an 
EHB ([Rosenberg et al ! l2004h . Finally, the very metal-poor remote halo cluster NGC 2419 also has 
an EHB. The one common property of these disparate clusters is that they are all among the most 
massive objects in the Galactic globular cluster system. 

Both metallicity spreads and He-enhanced populations imply multiple episodes of star formation 
and self-enrichment in some globular clusters, likely involving Type II supernovae and/or winds 
from massive AGB stars. This picture is consistent with the observed correlation with mass, as 
only the most massive clusters are likely to be able t o reta in significant amounts of ejected gas. For 



example, the scenario proposed bv lD'Antona et alJ ([20051 ) to explain the NGC 2808 main sequence 



and EHB proposes three distinct episodes of star formation spread over several hundred Myr: the 
initial burst at big-bang He abundance (Y ~ 0.24), followed by a second generation with Y ~ 0.4 
born from the winds of massive (~ 6 — 7M & ) first generation AGB stars, and later a third generation 
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with Y ~ 0.26 — 0.29 born from the winds of less massive (~ 3.5 — 4.5M ) AGB stars. Similar 
scenarios have been proposed for oj Cen, although these often also include the possibility that this 
cluster was formerly the nucleus of a now-defunct dwarf galaxy. Numerous problems remain to be 
solved with self-enrichment models, involving, for example, the retention and mixing of ejected gas 
within a cluster, the required first-generation IMFs (to get enough AGB stars), and heavy-element 
pollution from Type II supernove. 

We propose the following MODEST question and task: 

• In the next one to three years: can material enrichment take place as a result of multiple 
mergers in very massive globular clusters? 

• In the next three years: development of more comprehensive models of self-enrichment (includ- 
ing hydrodynamic modelling of accretion from interstellar gas within clusters) and multiple 
star- formation episodes within globular clusters. 

3.4 Observing stars escaping from globular clusters 

Measuring the properties of stars which are escaping from a particular globular cluster has the 
potential to tell us much about the internal processes in that cluster (eg, Gunn & Griffin, 1979). 
In particular, it would be of interest to try and find stars which have been ejected from the cluster 
with some significant velocity, as these offer a means of probing strong interactions between cluster 
members, such as three-body and four-body encounters (eg, Meylan et al, 1991). One possible way 
to locate such stars would be by measuring proper motions in a nearby target (e.g. M4) near the 
tidal radius or Lagrangian points. Once suitable candidates have been located, radial velocities 
could also be obtained. Because proper motion measurements require multi-epoch imaging over a 
significant baseline, this is necessarily a problem to be tackled on a time-scale of at least several 
years. Modelling can be utilized to predict the number of expected detections, and will be vital in 
attempting to constrain the processes which could produce the observed properties of any strong 
candidate high- velocity escapers. It is worth noting that for several Galactic globulars (e.g., Palomar 
5, Odenkirchen et al, 2003; Leon et al, 2000), large numbers of escaped stars have been observed 
in the form of tidal tails. These are members which have drifted through the Lagrangian regions 
with small relative velocity, and are more useful for probing the cluster's orbit about the Galaxy 
as well as the Galactic mass distribution. 

We propose the following MODEST questions: 

• In the next three years: can we locate new tidal tails belonging to any Galactic globular 
clusters? Searches utilizing deep wide-field imaging are presently in progress. 

• In the next three years: from a modelling aspect, can we predict the number of expected 
detections of fast escapers from a given globular cluster? Can we also constrain the processes 
which produce them, and predict their observational properties (e.g., the velocity distribution 
of fast escapers)? 

• In the next ten years: can we locate high proper motion stars near the tidal boundaries of 
nearby globular clusters? If so, radial velocities of any candidate stars need to be obtained. 
What can the observed properties of such objects tell us about the internal processes in their 
parent clusters? 



7 



3.5 Observational constraints on binary star populations in globular clusters 

Relatively little is known about binary star populations in globular clusters. Photometrically, we 
can begin to deduce probable binary candidates from the binary main sequence (eg, Rubenstein 
and Bailyn 1997). This is a region of the CMD lying above and to the red of the single-star main 
sequence. An unresolved binary consisting of two main sequence stars will have a combined colour 
somewhere in between the colours of the two components, and a magnitude brighter than that 
of the single-star main sequence at this combined colour. In principle, it is therefore possible to 
determine cluster binary fractions by observing the binary second sequence; however in practice 
this process is complicated by photometric errors, which mimic the main sequence spread due to 
binary stars, as well as crowding and field star contamination. Hence, to date binary fractions have 
only been measured at low significance in Galactic globular clusters via this method. Photometric 
variability surveys are also sensitive to some types of cluster binaries - in particular those which 
have the correct inclination to be eclipsing objects, plus those which are active in some manner, 
for example contact binaries (eg, Kaluzy et al. 1999). However, the relationship (if any exists) 
between active binaries and the global population is not known, so it is difficult, if not impossible, 
to infer properties of the normal binary population from the active one. 

Cluster binaries can also be detected via spectroscopic observations (eg, Pryor et al 1989 ; Hut 
et al. 1992 ). The idea is to obtain radial velocity measurements at multiple epochs. Binary stars 
should show large variations in the measured radial velocity due to the orbital motions of the two 
components. Such measurements have the advantage that they can provide information about the 
orbital period, providing the sampling is suffiently frequent. The disadvantage of this technique 
is that it requires repeated time-consuming observations, and is rather inefficient considering the 
binary fraction is expected to be of order 10% or less. 

A third way of probing a cluster's binary star population is by means of X-ray observations. 
With the launch of Chandra and XMM-Newton, many new X-ray sources have been detected in 
Galactic globular clusters. From all ROSAT observations, 57 X-ray sources were discovered in 
as many globular clusters (Verbunt 2001), but in 47Tuc Chandra already found over 300 sources 
(Grindlay et al. 2001, Heinke et al. 2005). From Chandra observations of 12 globular clusters, 
Pooley et al. (2003) found that the number of X-ray sources with an X-ray luminosity above 
4 x 10 30 ergs _1 scales with the collision number. This scaling was interpreted as evidence that 
these X-ray sources, which are expected to be primarily cataclysmic variables, are formed through 
dynamical interactions. From the X-ray colours and luminosities it is possible to identify and 
classify the sources containing an accreting neutron star (a low-mass X-ray binary), but for the 
large part of the X-ray sources optical identifications are necessary to discriminate cataclysmic 
variables from magnetically active binaries (RSCVn, BYDra, WUMa systems). Sofar, this has 
been done for 4 globular clusters; NGC 6752 (Pooley et al. 2002), 47Tuc (Edmonds et al. 2003ab), 
NGC6397 (Grindlay et al. 2001) and M4 (Bassa et al. 2004). These observations suggest that the 
majority of the X-ray sources with X-ray luminosities above a few times 10 30 erg s" 1 are cataclysmic 
variables, while the fainter sources are active binaries. From a comparison of these identifications 
in M4, 47 Tuc and NGC 6397, it is found that the number of such active binaries appears to scale 
with the (core) mass of the cluster, instead of the collision number (Bassa et al. 2004). Hence it 
seems likely that these systems evolved from primordial systems. 

We propose the following MODEST questions and tasks: 

• In the next year: can we make any robust estimates of binary star fractions in globular 
clusters for which high-quality main sequence photometry already exists? This will require 
sophisticated statistical treatment of observational errors to disentangle the binary second 
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sequence. Does the spatial distribution of binaries look like what is predicted from dynamical 
models of star clusters (eg, Ivanova et al 2005)? 

• In the next three years: the establishment of a detailed observing program of one nearby clus- 
ter (possibly M4), with the aim of determining the binary fraction, along with the distribution 
of mass ratios, periods, and separations for these objects. This will require both photometry 
and multiple-epoch spectroscopic measurements. In addition, can we start to establish a link 
between the active binary population in this cluster, and the overall population? This would 
allow us to start inferring the properties of binary populations in other clusters based on 
already-existing observations of the active systems. 

• In the next ten years: if the M4-type measurements can be extended to several other clusters, 
we can start to build a global picture of binary populations in globulars, including intra-cluster 
variations. How do the characteristics of the binary populations vary with cluster properties 
(e.g., mass, concentration, etc)? What can this tell us about internal cluster dynamics? 



3.6 Constraining initial conditions for globular cluster simulations 

In recent years, star cluster simulations have reached new levels of power and sophistication. For 
example, we are now in a position to run direct collisional iV-body models of objects at the lower 
end of the globular cluster mass function, incorporating much realistic physics (such as stellar and 
binary star evolution), for example see Hurley et al. 2005, where the open cluster M67 was modelled. 
With this type of modelling comes new challenges. One of the chief among these is the question of 
what initial conditions should be used for direct, realistic globular cluster simulations. Of course, 
the initial conditions adopted for any given run depend strongly on the system being modelled 
and the aims of the simulation. Nevertheless, it is important to develop a global understanding 
of how initial conditions affect subsequent evolution, in order that the most suitable starting point 
can be selected for any given simulation: can we constrain which initial conditions strongly affect 
subsequent long-term cluster evolution, and which are essentially irrelevant; or observationally, do 
we see any objects which will evolve into globular clusters like those in the Galaxy over the next 
Hubble time? Examples of modelling of globular clusters can be found in Phinney (1993), Druckier 
(1995), and Giersz & Heggie (2003). 

There are two ways to address this problem. From a modelling point of view, as more and more 
large-scale simulations are calculated, with varying initial conditions, it should become clear which 
parts of parameter space (especially covering the IMF and initial spatial and velocity structures 
such as mass segregation) strongly influence cluster evolution, and how. Observationally, detailed 
measurements of very young star clusters can help provide constraints on realistic initial conditions. 
Such measurements have alre ady been utilized in a number of studies - for example, the direct 
modelling of LMC clusters by ( Mackev et al. . 20061 ) involved initial conditions strongly constrained 



by the observed properties of young massive LMC objects such as R136, NGC 1805, and NGC 1818. 
We propose the following MODEST questions: 

• In the next year: begin to run simulations designed specifically to investigate the influence 
of initial conditions on the early, intermediate, and late-time evolution of globular clusters. 
Can we identify which are the most important initial conditions and which have little or no 
effect? Can we infer initial conditions from already-existing observations of young clusters? 

• In the next three to ten years: Can one run direct or near-direct simulations with the aim of 
demonstrating whether the super star clusters observed in starburst and interacting galaxies 
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are really globular cluster progenitors? If they are not, then what did the Galactic globular 
clusters look like initially? 



3.7 Are there black holes in globular clusters? 

Low-mass black holes (LMBHs) are expected as the end products of the evolution of stars pop- 
ulating the uppermost end of the IMF - i.e., those with M > 20M Q . If such remnants are 
formed without significant initial velocity kicks, the retention fraction is ex pected to be high and 



tormed without significant initial velocity kicks, tne retention traction is ex pected to be high and 
the holes should constitute a d ynamically important cluster sub-population (jKulkarni et al I ll 99,4 



Sigurdsson fc Hernauist . 1993h . Depending on the shape of the upper IMF, most if not all globular 



clusters are expected to possess this population of up to several hundred LMBHs early in their 
life. Within ~ 1 Gyr of formation, most LMBHs in a cluster have settled via mass segregation 
to form a centrally concentrated core. Eventually this core is sufficiently dense that multiple-hole 
interactions occur, resulting in the formation of BH-BH binaries and the ejection of single holes. 
Subsequent interactions harden the BH-BH binaries until eventually the recoil velocity is high 
enough for ejection. Several single holes are also expected to be ejected during this hardening pro- 
cess. Hence, it is thought that the LMBH population in a cluster completely depletes itself within 
a few Gyr. Nonetheless, the LMBH population is expected to inject significant amounts of energy 
into the stellar core in a cluster before depletion, both through the dynamical scattering of stars 
and the removal of BH mass from the cluster centre. A-body simulations show that this influence 
is in many cases enough to significantly alter (expand) the structure of the stellar core. Therefore, 
LMBHs likely represent an important (and o ften neglected ) dyna mical influence in the early and 
intermediate phases of star cluster evolution ( Mackev et all l2006h . 



Intermediate-mass black holes (IMBHs) in clusters are interesting because it is thought that 
such objects may represent the seeds of the super-massive black holes (SMBHs; M > 1O 6 M ) which 
are inferred to exist both in high-redshift galaxies (where they are believed responsible for quasars 
and AGN), and in the local universe at the centres of our Galaxy and M31. Stellar dynamical 
simulations suggest that IMBHs can be formed in very dense young globular clusters, via the 
process of runaway merging. In such clusters, the core-collapse timescale for the most massive stars 
can be shorter than their main sequence lifetimes. Core collapse may initiate a rapid sequence of 
direct collisions between stars, le ading to the production of a merged o bject (possibly an IMBH) of 
mass ~ 0.1% of the cluster mass ( Portegies Zwart &; McMillanl . l2002bh . It is likely the inspiral and 



destruction of a suitable cluster near a galactic centre (e.g., the Arches or Quintuplet in our Galaxy) 
may seed or contribute to the growth of super- massive black holes (M > 10 6 Mq) as described in 
Section 4. The possibility of detecting IMBHs in globular clusters is thus both highly relevant and 
intriguing. 

To date, we possess only indirect and/or debated evidence for IMBHs or populations of LMBHs 
in globular clusters. The presence of a ~ 2500Mq IMBH in the nearby globular cluster M15 
dGerssen et alll2002h and a 1.8xl0 4 M o IMBH in the massive stellar cluster Gl in M31 (|Gebhardt et al 



200.4) have been i nferred from HST measurements. However, these detections are contested by 



Baumgardt et al.l (|2003al lbD. whose A-body modelling suggests that neither detection requires the 
presence of an IMBH - that is, each set of measurements can seemingly be explained by models 
with large central populations of stellar remnants such as neutron stars and white dwarfs. It has 
been suggested that detected X-ray emission from Gl may be the result of accretion of gas by a 
central black hole (Pooley & Rappaport, 2006). The more general question of which globular clus- 
ters may contain IMBHs has also been considered in Baumgar dt et al (2005). The only e vidence 
for LMBH populations in globular clusters is the observation ( Mackev fc Gilmorel . 2003al lbh that 



intermediate age clusters in the Magellanic Clouds possess a wide range of core sizes, and that this 
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range apparently correlates with cluster ag e. This trend can, at least in part, be explained by the 



dynamical influence of LMBH populations ( Mackev et all l200fil ) . 



The main problem in detecting the presence of IMBHs (or centrally concentrated LMBH pop- 
ulations) in globular clusters is that these objects have only small spheres of dynamical influence: 
radius ~ 0.1 pc for a 2000Mq IMBH in a typical cluster core. At present therefore, measurements 
are limited by resolution. This problem will likely be solved on a time-scale of at least ten years, 
with the advent of 50 — lOOm-class telescopes and functional adaptive optics at visible wavelengths. 
Togethe r these woul d permit full (3D) dynamical studies of globular clusters out to the Magellanic 
Clouds (|HoolJ . l2005h . allowing secure detections of IMBHs in globular cluster cores, as well as the 



possibility of investigating the effects of LMBH populations in intermediate-age clusters. In the 
meantime, ever more realistic modelling can help place constraints on both IMBH and LMBH for- 
mation in clusters, as well as what observational signatures should be expected from the presence 
of such objects. 

We therefore propose the following MODEST questions: 

• In the next year: using iV-body modelling, can we make any new predictions about the 
observational signature(s) of a population of LMBHs in a globular cluster? 

• In the next three years: can we resolve the disagreement between the observations of objects 
such as M15 and M31-G1 which infer the presence of IMBHs, and iV-body modelling which 
suggests IMBHs are not required in order to explain the observed dynamics? Can we refine 
models of stellar evolution to make concrete predictions about whether we should really be 
expecting IMBHs or large populations of LMBHs in globular clusters? 



4 Galactic nuclei 



4.1 The final parsec problem 



It is generally accepted that hierarchi cal models best expla i n the formation of structures in the Uni 
verse , down to the size of a galaxy ( White fc Reesl . 197Sl : Kaufimann et al. . 1199.4 Springel et al 
200.51 ). This means that in their lifetimes galaxies typically merge with one or more other galax- 



ies. A good example of this is the A ntennae "galaxy", which actual ly consists of two colliding 



galaxies, NGC 4038 and NGC 4039 dWhitmore fc Schweizerl. Il995bh. Almost a ll galaxies ap 



^ar to have a central supermassive hWk hoi. fSMBH. krrarese Ford id for" a re v^Y 



hence mergers of galaxies can eventually lead to me rgers of the SMBHs (e.s.. iMenou et all l200ll : 
Haehnelt fc Kauffmannl 120021 : IVolonteri et alll2003h . At the beginning of the evolution, dynamical 
friction makes the orbits of the SMBH decay, so that they sink down to the centre of the merging 
system. Strong interactions with surrounding stars coming from the stellar system in which the 
SMBHs are embedded remove energy and angular momentum from the SMBHs after they have 
formed a bound binary system. These stars are re-ejected into the stellar system with an increased 
kinetic energy and thus the semi- major axis of the SMBH binary shrinks. The rate of shrink- 
ing slows down after the SMBHs are close enough that they are more massive than the enclosed 
stellar mass. This typically happens at a separation ~ 0.1 — 1 pc, significantly larger than the 
0.001 — 0.01 pc needed so that gravitational radiation alone can cause the binary to merge within 
a Hubble time. The "final parsec problem" ( Begelman et al. . 198fll : Milosavljevic fc Merrittl . 20031 : 
Yul l2002l l thus consists of identifying processes that can bring the binary separation from ~1 pc to 
the realm of significant gravitational radiation. The efficiency of such processes has major implica- 
tions for the growth and mergers of SMBHs, galaxy evolution, and sources for future space-based 
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gravitational wave detectors such as the Laser Interferometer Space Antenna (LISA). For a com- 
prehensive review of th e vario us aspects of formation and evolution of binary MBHs we refer to 
Merritt & Milosavlievio! (|2005h . 



In a galactic nucleus that has negligible amounts of gas and that has an axially symmetric 
potential, individual stars in the nucleus essentially conserve their angular momentum from orbit 
to orbit. Thus, after the stars whose pericentres take them near the SMBH bin ary are ejected 
from the system via a gravitational slingshot (hence the "loss cone" is emptied; see Frank Sz Reed 
further tightening of the binary requires that distant two-body encounters between stars send 
some of them on orbits radial enough to interact with the binary. The timescale for this two-body 
relaxation can be billions of years, but recent work based on the empirica l relation between the mass 
of an SMBH and the velocity dispersion or mass of i ts host galaxy bulge dFerrarese Merrittl . Eoool: 
Gebhardt et~ail l2000l : iMerritt fc Ferraresel . l2001al lbl: iTremaine etail 120021 ) suggests that SMBH 
binaries with total mass Mrrh < 10 7 M Q can be hardened to merger in less than a Hubble time 



( Milosavlievic Sz Merrittl . I2003T ) . This is the realm relevant to low- frequency gravitational wave 
detectors su ch as LISA. In contrast, it is still unclear whether a combination of other factors such 
as gas drag dEscala et al.l.l200,4lMilosavlievic fc PhinnevL l2005h o r triaxiality of the galac tic nuclear 
potential (jHollev-Bockelmann et al 1 120021 : IPoon &r. Merrit4l2002l : IMerritt Poorl l2004h suffices to 
produce efficient mergers for more massive SMBH binaries. 

Ideally, one would like to ensure inclusion of all relevant physical effects with a direct-summation 
N-body treatment of galactic centre dynamics. However, the actual number of stars in the central 
few parsecs of a galaxy is 10 7-8 , which is too much to simulate in this way even for special- 
purpose supercomputers such as the GRAPE-6, which reach a computational power of 64 Tflops 1 . 
Using fewer particles decreases the effective two-body relaxation time and therefore introduces an 
artificially high rate of binary hardening. 

A possible solution is parallel usage of a direct-summation N-body code on a cluster of special- 
purpose GRAPE-6 nodes. However at the present time, there is no N-body code which treats close 
encounters rigourously (through Kustaanheimo-Stiefel two-body as well as chain regular isation) 
and is both adapted for the use on GRAPE-6 and fully parallelized (jAa.rsethl . IT99I l2003ah . 

Another possibility is that the real dynamics are not so sensitive to the number of stars. For 
example, it has been proposed that non-axisym metries could help the binary shrink much faster 
thanks to the chaotic nature of the stellar orbits ( Hollev-Bockelmann et al. , 20021 : Poon &: Merrittl . 
I2OO2I : IMerritt fc Poonl . 120041 ) . Recent work based on this approa ch suggests that the rate of orbital 
decay is roughly independ ent of the total number of particles ( Hollev-Bockelmann fc SigurdssonL 
20061 : iBerczik et all 120061 ). However, caution is appropriate because extrapolation to the much 



larger number of stars in real galaxies requires theoretical scalings which are not fully understood. 
N-body simulations should be therefore envisaged as a source of encouragement and motivation 
rather than solid and robust proofs. 

For the MODEST questions on this subject we propose: 

• In the next year: observationally, what is the dynamical state of a galactic centre just after 
In particular, what is the rotational structure, and how triaxial are the 



? 



a major merger 
centres? 

In the next three years: what is th e influence of nuclear rotation on the dynamics of a SMBH 



binary? Initial investigations, e.g.. IBerczik et alJ (|2006h . suggest that the influence could be 



substantial, but this needs to be coupled with observations. 



1 http : //grape . astron. s .u-tokyo . ac . jp/grape/ 
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In the next ten years: what is the evolution of the eccentricity of an SMBH binary from its for- 
mation until it becomes detectable with low-frequency gravitational radiation detectors such 
as LISA? This will req uire incorporation of ro t ation, triaxiality, and the effects of resonances 
on dynamical friction (|Tremaine fc Weinber and substantial input from analytical 

treatments, direct N-body summation techniques, and more approximate approaches. The 
potential payoff is that any residual eccentricity detected with LISA or similar instruments 
might then be used to untangle important elements of the mergers. 



4.2 Stellar dynamics around a massive black hole 

The central SMBH and the stellar system interact through many channels in addition to the smooth 
gravitational potential. For example, stars can produce gas to be accreted on to the SMBH, through 
normal stellar evolution, collisions, or disruptions of stars by the strong central tidal field. These 
processes may contribute significantly to the mass of the SMBH. Tidal disruptions trigger phases of 
bright accretion that may reveal the presence of an SMBH in an otherwise quiescent, possibly very 
distant, galaxy. Collisions may create observationally peculiar stellar populations. Also, stars too 
compact to be tidally disrupted are swallowed whole if they are kicked directly through the horizon 
("direct plunges") or progressively inspiral down to a relativistic unstable orbit through emission 
of gravitational waves (GWs). The latter process, known as an "Extreme Mass Ratio Inspiral" 
(EMRI) will be one of the main targets of LISA. 

Many different numerical schemes have been applied to the simulation of galactic nuclei hosting 
an SMBH. Most of them rely on the assumptions of an isolate d, spherical system in dynamical 
equilibrium (e.g., direct integration of Fok ker-Planck equation by Murphy et al. 199ll . Monte-Carlo 
methods by Duncan fc Shapirol 19821 and Freitag fc Bend 20021 and gas-dynamical treatment by 
Amaro-Seoane et alJ 20041 ). Inthese cases, only "collisional" effects can bring stars on to the loss- 



cone, i.e., the very elongated orbits which allow close interaction between a star and the SMBH. 
These effects include (diffusive) 2-body relaxation, large-angle scatterings, direct collisions, and 
resonant relaxation. The approximate methods just mentioned generally include only diffusive 
relaxation and, in rare cases, collisions and large-angle sca tterings. There have only recently been 
direc t N-body simulations of clusters with a central object ( Baumgardt et al. . 2004al bl: Preto et al. . 
120041 . in particular). 

As with binary SMBH simulations, for these applications N-body simulations are invaluable 
because they dispense with the necessity of most approximations but it seems unlikely that they 
will completely supersede other approaches in the next few years. Their most obvious limitation 
is the steep computing time scaling icPU °^ ^part w hich currently limits the number of particles 
to about 10 6 , short of the 10 7 - 10 8 stars in even small galactic nuclei. More fundamentally, 
the community is missing an algorithm suitable to tackle this particular class of problems. The 
integration of millions of nearly Keplerian orbits with the Hermite scheme used in usual N-body 
codes causes spurious change in the orbital constants and may, for instance, lead to incorrect star- 
MBH interaction rates due to eccentricity increase. Therefore, numerical errors dominate the effects 
of the actual star-star perturbations. 

The development of an N-body regularisation scheme to integrate the motion and mutual pertur- 
bations of a large number of light objects orbiting the same massive body has become a challenging 
priority. The chain scheme is not suitable here because it do es not allow reg ularisation of the 

for technical back- 



priority. Ine cnam scheme is not suitable nere because it do es not allow reg 
interaction of the SMBH with each close star simultaneously (see Aarsethl 2003b f 



ground). Some inspi r ation may be provided by t he symplectic codes in use in planetary dynamics 
( Wisdom Sz Holmanl . 199ll : I Wisdom et al. , 19961 ) . These codes are optimised for nearly-Keplerian 
orbits with a single dominant object, and bound the energy error over large numbers of orbits. 
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They may therefore have applicability to the central regions of galactic nuclei that are dominated 
by an SMBH. 

It is, however, likely that new development will be necessary for dynamics around an SMBH 
because excellent accuracy in all the orbital elements, not just the semimajor axis, may well be 
necessary to preserve the proper interactions with narrow resonances or to include the correct effects 
of processes such as general relativistic precession. 

The inclusion of relativistic contributions t o dynamics is especially im portant for the proper 
treatment of effects such as resonant relax ation 
19621 : iLidov fc ZiglirJ Il97fil : llnnanen et al 



1997 



Rauch fc Tremainel 199(f) or Kozai cycles ( Kozail . 
Miller fc Hamilton! I2OO2I ) . which depend on the 



persistence of certain phase relations over hundreds of orbits or more. Even relatively small effects 
can have an influence over this many orbits, hence precision of integration and inclusion of relativity 
are at a premium. In addition, estimates of the distributions of mass ratio, eccentricity, etc., are 
needed to construct reliable template banks for detection of the gravitational radiation from EMRIs. 
Simulations based on purely Newtonian schemes may well lead to completely incorrect results 
for these purposes. General relativistic effects up to the 2.5 post-Newtonian order (and hence 
including both preces sion and radiatio n reaction) have b een implemented in the codes HNBody 



( Giiltekin et al. , 2006) and N-body6++ (|Kupi et al.l . l2006l ). but the applications are in their infancy. 



Whatever the eventual solution, it is likely that such N-body codes will initially lack the effi- 
ciency to treat more than 10 4 — 10 5 stars, corresponding to the inner part of the influence region 
of a realistic nucleus. An important first step would be to consider the central SMBH as a fixed 
particle, i.e., treat it as an external potential. Neglecting the random motion of the SMBH may 
have consequences for, e.g., the rate of tidal disruptions because this motion should allow loss-cone 
replenishment. However, this should only become an issue when a region larger than the influence 
radius can be simulated. Furthermore, this idealisation will prove useful for comparison with ana- 
lytical and approximate numerical approaches which generally rely on it. It is also the best way to 
assess later the role of the SMBH motion by comparison with more realistic simulations. 

To pave the way for N-body studies and to complement them, more approximate but much 
faster and more flexible methods are invaluable. Monte Carlo (MC) statistical approaches seem 
ideal because, being based on particles, they make it easier to follow the evolution of individual 
orbits and to include individual star-star or star-SMBH interactions. A prom ising avenue would be 
to combin e aspe cts of the MC approach of lHopman Sz Alexandei (1200.4 2006) to that first pioneered 



by Henonl ( 197.4 ) for globular cluster dynamics. iHopman fc Alexanderl followed the orbital evolution 
of test particles due to diffusive and resonant relaxation and to GW emission in the Keplerian 
potential of the SMBH assuming a fixed stellar background. The Henon approach would evolve 
the stellar distribution self-consistently -thus obtaining the correct m ass-segregation effects- but 
its only application to galactic nuclei so far ([Freitag &: Bend . l2002h lacks the ability to resolve 
the dynamics of stars in or near the loss cone on sat isfactorily short timescales. Shapiro and 
collaborators ( Duncan &: Shapircl . Il98l IShapiroL Il985l) developed an MC code which represents 
the cluster as a set of spherical shells like in the Henon scheme but where the effects of relaxation is 
computed by explicit integration of diffusion coefficients like in direct Fokker-Planck codes, rather 
than pairwise interactions between particles. This method had the important advantages of allowing 
one to increase the resolution in the central regions and to follow particles very close to the SMBH 
on orbital timescales, but it has only been used for single-mass cases. Such an algorithm, if it can 
be extended to a stellar mass spectrum, would be an extremely useful tool. 

Theoretical predictions for ra tes and characteristics of EMRI events have proven quite p roblem- 

for a quick review and Hopman fc Alexander 200.4 200fil for new 



atic so far (see| 

developments). Even in the standard case where only diffusive relaxation on single stars is con- 
sidered, different authors find results scattered over such a large range that it is not clear whether 
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only a few events will be detected with LISA ( Hopman fc Alexander , 2005h or whether they will 
turn out to be an embarrassment of ri ches, preventing the individual detection of each other and of 
other sources ( Barack fc Cutler . 20041 ). Recently, non-standard processes to bring compact objects 
very close to an SMBH have been sug gested. These include the accu mulation of r ed-gia nt cores 



by tidal peeling ( Davies &: King . 2005 ) , tidal separation of b inaries ( Miller et al. . 2005h . stellar 



formation in an accretion disk ( Levin 



200,4 iNavakshinl . 1200.4 ) and dissipative interactions with a 



disk (S ubr et al.l . 120041 ) . In all these cases, the orbital evolution will start being dominated by GW 



emission (as opposed to relaxation) at a smaller semimajor axis and much smaller eccentricity than 
in the "standard" case and the EMRI should be very nearly circular in the LISA band. The study 
of these various types of EMRIs will be one of the main applications of the numerical methods 
envisioned here. 

Noticeably each of these processes was (also or uniquely) offered as a way to explain the origin of 



the young massive "S" stars orbiting SgrA* (e.g.. iGenzel et al. 1. 1200,4 . Hence it is no coincidence 
that, unless they somehow get natal kicks comparable to their orbital speeds of thousands of 
kilometers per second, the remnants of the "S" stars will have the appropriate orbital parameters 
to become EMRIs, a fact with strong bearing on LISA detection rates if the situation around Sgr A* 
is typical. 

It has also been suggested that intermediate-mass black holes (IMBH, with masses ~ 10 2-4 -M©) 
formed in young massive star clus ters within ~ 100 pc of the centre of a galaxy could sink to the 
centre and merge with an SMBH (|Millerl 1200.4 IPortegies Zwart et all 120051 : iMatsubavashi et all 
20051 ) . If this happens it would be an extremly stro ng sou r ce for LISA, with unique potential for 
mapping the spacetime around a rotating SMBH ( Miller . bo05h . Currently, however, there are 
many uncertainties about the various steps in the sinking process, from the settling of a cluster in 
a galactic nucleus to the stripping of that cluster to the processes that allow an IMBH to merge 
with an SMBH (for example, is there any stalling and if so, will other IMBHs come in and cause 
mutual ejection?). 

Finally, while structures such as triaxial bulges, bars or stellar discs are common on scales 
of 100-1000 pc, the influence of non-sphericity at small and intermediate scales on the structure 
and evolution of the nucleus has been little explored. The existence of a large fraction of "cen- 
trophilic" (box and chaotic) orbits in triaxial st ructures has the potential of boosting the rate of 
star- S MBH interactions by orders of magnitude dHollev-Bockelmann et ah . 20021 : Poon &: Merrittl . 
20021 : iMerritt PoonL l2004 iHollev-Bockelmann SismrdssonL |200fih . For EMRIs, though, it is 
not clear whether such orbits, with very large initial semimajor axis and eccentricity, have a chance 
to shrink to L/5^4-detectable frequencies without being perturbed into a direct plunge or a wider 
orbit. 

We therefore propose the following MODEST questions: 

• In the next year: are the basic codes used to calculate EMRI rates consistent with each 
other? As proposed to us by Richard Mushotzky, we suggest that a precisely defined test case, 
accessible to direct N-body summation methods as well as to statistical approaches, should be 
simulated by several independent groups. Such comparison has proved very enlightening in 
the case of th e "collaborative experiment" in cluster dynamics organised by Douglas Heggie 
dHeggiel I2003I . see also http : / /www . manvb ody . org/mo d est/W G/ wg7 . html) or in the field of 



cosmological hydrodynamics simulations ( Frenk et al. . 19991 ). For example, current codes 



could treat a cluster of 10 4 point-mass single stars hosting a central SMBH with a mass 
1% of the total. Several quantities could be compared, including the time evolution of the 
capture rate. The observed similarities and differences could help guide further treatments. 
For the number of EMRIs to be significant with such a particle number, the cluster needs to 
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be made more compact than any known real system, to boost relativistic effects relative to 
relaxation. Although lacking physical realism, this setting will allow one to test and calibrate 
the approximate methods by comparison with direct N-body. 

In the next three years: what are the capture and tidal destruction rates implied by the actual 
distribution of stars around the SMBH in our own Galaxy? The radial dependence of number 
density is reasonably well constrained within ~0.01 pc of the centre, a region containing some 
10 4 ~ 5 stars, so a specific simulation over the ~ 10 9 yr relaxation time would be informative. 

In the next ten years: what is the true influence of nonaxisymmetry at large distances on the 
inner few parsecs, where EMRIs interact? Does this lead to large rates of LISA detections, or 
does it instead produce direct plunges? Before direct N-body methods are able to deal with 
> 10 7 particles, the relaxational dynamics of non-spherical systems could be studied with 
hybrid schemes borrowing from "collisionless" N-body and Fokker-Planck o r MC codes, an 
option still virtually unexplored (with the exception of Johnston et al. 19991 ) . As a separate 



but related matter, what are the processes that lead to an IMBH-SMBH merger, and what 
is the expected rate of LISA detections? 

4.3 Stellar collisions 

Galactic nuclei are one of the few environments in which collisions involving single stars should 
occur on a relatively short timescale. For instance, within ~ 0.03 pc of Sgr A*, a 1M Q main 
sequence (MS) star should experience, on average, one collision in less than 10 Gyr. For a giant, 
this timescale is reduced to a few 10 7 yr. Although collisions probably do not strongly influence 
the stellar dynamics, they are of great interest as a way to produce unusual stellar populations. 
They have been suggested as the cause of the apparent paucity of giants in the vicinity of Sgr A*, 
although giants irradiated by the X-ray radiation of Sgr A* may actually masquerade as massive 
MS stars ( Jimenez et al. . 2006^ . 



Collisions in galactic nuclei occur at relative speeds of a few 100 km s" 1 or higher, making 
mergers an unlikely outcome for low-mass stars. In particular the possibility of growing "super 
blue stragglers" through a sequence of collisions seems excluded. The m ass and energy loss for such 
high- velocity collisions between MS stars has been studied exhaustively ( Freitag fe Benz . 20051 ) but 



much remains to be done for the more likely case of a collision between a giant and a more compact 
object. In relatively small galactic nuclei (typically hosting an SMBH less massive than 1O 7 M ), 
collisions involving stellar BHs are also of special im portance because m ass segregation probably 
concentrates these objects around the central SMBH dFreitag et al l bnod ). 



High-velocity collis ions between a giant and a smaller star were computed, using SPH, by 



xy coins 

Bailev fc Daviea (|l999h who found that during a typical collision the impactor, flying through the 



giant's envelope, causes only relatively little mass loss; the giant is likely to recover on a short 
(thermal) timescale. Complete removal of the envelope can only happen if the smaller star is 
captured and a common-envelope (CE) system is formed, an outcome too rare t o explain the 



dearth of giants at the Galactic centre. On the other hand, iDavies et al.l (|1998I ) showed that 
collisions between giants and binary stars may be more efficient at depleting the giant population, 
either by creating CE systems or by ejecting the giant's core from its envelope, if binaries are 
commo n enough. Howeve r, except for the detection of transient X-ray sources at < 1 pc from 
Sgr A* ll.Muno .1 all 1200.4 ). very little is known about binary populations in galactic nuclei. 



Our MODEST questions are: 
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In the next year: what are the dynamics of binaries in galactic nuclei? The central goal is 
to study the survival of binaries in an environment with a much higher velocity dispersion 
than exists in globular clusters. It will be particularly important to study this question using 
stellar dynamical simulations with a large number of particles (see Sec. 14.2)1 . 

In the next three years: what is the evolution and appearance of a giant star whose envelope 
has been partially removed by a collision (or a strong tidal interaction with the SMBH, 
see Di Stefano et al. 2001 : Davies fc Kinj 20051 ) . Also, what is the evolution of common 



envelope binaries formed through red giant collisions? These systems may be the progenitors 
of compact binaries, possibly explaining (some of) the X-ray sources observed around Sgr A*. 

In the next ten years: what are the hydrodynamic and possibly magnetohydrodynamic results 
of collisions between giants and smaller stars, and between compact objects (especially stellar- 
mass black holes) and extended stars? For the former, simulations need to cover a much 
more extended region of parameter space (masses, evolutionary stage of the giant, relative 
velocity and impact parameter) than published so far. For the latter, it will be important 
to understand how damaging the collisions are, and how much mass the compact star can 
accrete. 
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